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String breaking with dynamical Wilson fermions 
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We present results of our ongoing determination of string breaking in full QCD with Nf = 2 Wilson fermions. 
Our investigation of the fission of the static quark-antiquark string into a static-light meson-antimeson system 
is based on dynamical configurations of size 24'^ x 40 produced by the TxL collaboration. Combining various 
optimization methods we determine the matrix elements of the two-by-two system with so far unprecedented 
accuracy. The all-to-all light quark propagators occurring in the transition element are computed from eigenmodes 
of the Hermitian Wilson-Dirac matrix complemented by stochastic estimates in the orthogonal subspace. We 
observe a clear signature for level-splitting between ground state and excited potential. Thus, for the first time, 
string breaking induced by sea quarks is observed in a simulation of 4-dimensional lattice-QCD. 



1. INTRODUCTION 

The fission of the flux string between two heavy 
quarks, is a key problem of lattice quantum chro- 
modynamics (QCD). While it is fundamental to 
a proper understanding of QCD, after nearly 25 
years of lattice QCD we are still lacking a con- 
vincing demonstration of string breaking through 
dynamical quarks Jj. 

Of course, extremely precise studies of the po- 
tential have been carried out, however, even at 
separations r = \R\a much larger than the antici- 
pated breaking distance of about 1 fm, still a lin- 
early rising potential is found. This phenomenon 
is attributed to the insufficient overlap of the Wil- 
son loop with the final state at large r. 

As first suggested back in 1992 by C. Michael, 
progress in the conceptual understanding of 
string breaking on the lattice could be achieved: 
in Ref. 0, QCD was mimicked by the the 
SU(2)-Higgs model using a two channel setting. 
The authors could confirm the interpretation of 
string-breaking as a dynamical level-crossing phe- 
nomenon between string-type and meson-type 



states. They observed a pronounced gap of the 
energy levels at the location of the crossing. 

A main hindrance of a straightforward appli- 
cation of this procedure to 4-dimensional QCD 
is the computation of all-to-all propagators to 
achieve a proper signal for the heavy-light con- 
tributions. The elements of the 2x2 mixing 
matrix — three of which are composites of Wilson 
lines and fermion propagators — must be aver- 
aged over all spatial lattice sites and dimensions. 
This requires the computation of the inverse of 
the Hermitian Wilson-Dirac matrix, Q^^. As the 
dimension of Q is about 0(10''), the full inverse 
cannot be represented even on high-end systems. 

In our investigation, we employ various "smear- 
ings" like APE smearing, Wuppertal smearing, 
fat temporal links and a new hopping parameter 
acceleration. The latter is shown to enhance the 
low eigenmode dominance. Furthermore, with 
the truncated eigenmode approximation (TEA) 
|3| , rendered exact by stochastic estimates on the 
orthogonal subspace, all-to-all propagators can be 
computed efficiently. 
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2. MODELLING THE TRANSITION 

Close to the location of QCD string breaking, 
we deal with a mixture of two states, the initial 
heavy QQ-pair at small separation r, represented 
by the standard Wilson-Loop W{R,T), and a fi- 
nal state at large separation r, which is given by 
two heavy-light mesons. The interesting physics, 
i.e. the fission of the flux-tube, takes place at in- 
termediate r. 

In the two-channel approach the correlator is 
represented as the following 2 x 2-matrix 
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Ci^i is the correlation function of the static QQ 
pair, namely the Wilson loop, 6*2,2 is the corre- 
lator of the heavy-light pair, BB, and Ci,2 the 
correlator of the transition matrix element. The 
light quark propagators in the generalized Wilson 
loops are represented by wavy lines. 

The exponential decay of the correlator is found 
by diagonalization of the real 2x2 transfer ma- 
trix, with the lower eigenvalue being the ground 
state energy of the system. 

As demonstrated in Ref. [2, characteristic sig- 
natures are required in order to "prove" string 
breaking: 

1. The energy levels show crossing with the 
asymptotic ground state reaching the same 
mass as the system of two separated heavy- 
light mesons. 

2. At the location of crossing, dynamical 
string breaking manifests itself in a pro- 
nounced energy level-splitting. 

3. The transition matrix element C12 will stay 
finite for a — *■ 0. 

3. SETTING 

The basis for our investigation of string break- 
ing are the Nf = 2 dynamical configura- 
tions produced by the SESAM and T^L groups. 
SESAM generated 8 ensembles of about 200 



decorrelated configurations each on 16'^ x 32 
lattices at two different /3-values, 5.6 (k = 
0.156,0.1565,0.157,0.1575,0.1580, 0.56 < ^ < 
0.83) and 5.5 (k = 0.158,0.159,0.1596,0.1600, 
0.67 < ^ < 0.85). These ensembles have been 
analyzed for string breaking in order to determine 
a (3 — K combination with small enough sea quark 
mass still being unaffected by finite-size effects. 
We have chosen the pair (3 = 5.6, k — 0.1575. 

With a Sommer scale of rg = 5. 892 (27) a at 
r = 0.5 fm the lattice spacing is a = 0.085 fm or 
fl-i = 2.32 GeV. ^ = 0.704(5) suggests a degen- 
erate quark mass close to the strange quark mass. 
For the 16'^ x 32 system, both the pion window of 
La ~ 4m^^ and the string breaking window with 
V2L/2rsB « .96 < 1 and V3L/2rB « 1.2 > 1 
are not large enough to comfortably accommo- 
date the string breaking experiment. We chose to 
use our T^L configurations generated on 24^ x 40 
lattices. Here the windows are La « 6'm~^ and 
V2L/2rsB « 1.2 > 1, V^L/2rB w 1.4 > 1. 
183 decorrelated configurations were analyzed for 
Wilson loops, 20 decorrelated configurations were 
analyzed for fermion lines. 

4. OVERLAP ENHANCEMENTS 

We use a set of mutually optimized methods to 
improve the ground state overlap which reduce 
the statistical error by a least a factor of 5: 

APE smearing. APE smearing is carried out 
on the gauge fields with a = 2.0 and N = 50: 



VSU(3) [U^{X) (2) 

^^(c/.(f)C/,(f + *)t//(f + .?) 
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Wuppertal smearing with (5 = 4 and linear 
combinations of Nw — 20, 40, 50: 
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Fat temporal links (HYP) a la Hasenfratz and 
Knechtli can lead to an exponential improvement 
of the signal quality. 
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Hopping parameter acceleration (HPA). 

The fcth power of the Wilson hopping term, 
(kD)'^, with M — 1 ~ kD, connects sites at dis- 
tance k and less. Multiplying = 1 + {kD) + 
{KDf + {KDf ... by [nDf leads to M-^inDf = 
M^^ — '^'^Ii (kDY ■ Thus, at distance fc, the first 
A; — 1 terms do not contribute, and the low eigen- 
mode dominance is enhanced, since with M\ip) = 

X\^P) we have M-^{KD)''\'tP) = 

Furthermore, we improve by the Trun- 
cated Eigenmode Approximation (TEA) 
[3]. TEA is made exact by means of stochastic 
estimates (SET) in the orthogonal subspace. In 
Fig. n we compare the standard SET- method (no 
mode computed via TEA) on the smeared gauge 
fields with the improvement through HPA and 
TEA (200 modes). 




Figure 1. Improvements by HPA and TEA. 

5. RESULTS 

Fig. 121 shows intermediate results of our ongo- 
ing evaluation of T^L configurations that have 
been obtained with a preliminary form of the fat 
link static action. Hence the asymptotic ground 
state is just slightly higher than twice the heavy- 
light mass. The error bars of the zoomed plot 
are smaller than the symbols. Level-splitting is 
achieved at the crossing point for the first time 
in 4-dimensional QCD with dynamical fermions, 
demonstrating string breaking. 
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Figure 2. Appearance of a level-gap. 
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